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Abstract 


Form factors of composite operators in the SL(2) sector of M = 4 SYM theory are 
studied up to two loops via the on-shell unitarity method. The non-compactness of 
this subsector implies the novel feature and technical challenge of an unlimited number 
of loop momenta in the integrand’s numerator. At one loop, we derive the full minimal 
form factor to all orders in the dimensional regularisation parameter. At two loops, we 
construct the complete integrand for composite operators with an arbitrary number of 
covariant derivatives, and we obtain the remainder functions as well as the dilatation 
operator for composite operators with up to three covariant derivatives. The remainder 
functions reveal curious patterns suggesting a hidden maximal uniform transcendental- 
ity for the full form factor. Finally, we speculate about an extension of these patterns 
to QCD. 
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1 Introduction and summary 


Recent years have seen increasing interest in the study of form factors. In particular, M = 4 
super Yang-Mills (SYM) theory in four spacetime dimensions furnishes a rich setting for 
studying their properties. Form factors in M = 4 SYM theory were initially investigated 
in 1985 [1]. Only recently, this subject was reconsidered, first at strong coupling 12, 3] and 
then at weak coupling [4, 5], and has subsequently aroused a lot more attention [6-36]. For 
recent reviews, see also the theses [37, 38]. 

An n-point form factor describes the interaction of n on-shell particles with momenta p; 
(p? = 0) and a gauge-invariant local composite operator O(x). Translated into momentum 
space, it takes the form 


Fon IT (1.1) 


where in this paper all form factors are considered to be colour-ordered super form factors 
and q is off-shell (q? £ 0). The motivation for studying form factors is multifold. Firstly, 
given the tremendous progress on on-shell techniques for scattering amplitudes (see e.g. [39- 
41]), it is desirable to transfer these powerful methods to the computation of more general 
quantities. Form factors combine external on-shell and off-shell objects and are thus ideally 


suited for this purpose. Showing a rich structure, they also provide new testing grounds 
for the further development of on-shell techniques. Secondly, forming a bridge between 
on-shell and off-shell observables, the study of form factors should yield novel insights 
into the hidden symmetries of the underlying theory. Moreover, we should emphasise that 
form factors are interesting in their own right. Sudakov form factors, for instance, have 
played a key role in the understanding of infrared singularities in gauge theories 142-451. 
Furthermore, they are closely related to Higgs production via gluon fusion as well as to 
Drell-Yan production in the standard model 146, 47]. 

The study of form factors of non-protected operators, as vvell as their connection to 
the spectral problem of planar W = 4 SYM theory, was recently pushed forward in 119— 
21, 28, 32, 36]. In particular, the form factor remainder was made a central concept and it 
was shown how to extract the dilatation operator from form factor data. While previous 
studies mostly focused on individual operators or compact sectors, in this paper we consider 
operators in the non-compact SL(2) sector. The crucial new feature of non-compactness 
brings technical challenges and reveals novel properties. 

Operators in the SL(2) sector of M = 4 SYM theory involve a single type of complex 
scalar X and an arbitrary number of insertions of covariant derivatives D" projected onto 
a lightlike vector T”: 


L : 
Dt)” 

on =i | J] oR x], D+ =D" r. (1.2) 
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The length of the operator, i.e. the number of scalars X, is denoted by L. The configuration 
of a certain state is thus completely specified by the vector 


A=(n1,...,nz), (1.3) 


which is identified with its images under cyclic permutations. We refer to the total number 
of covariant derivatives, M = —: nj, as the total spin or the total magnon number. In 
what follows, we will also speak of form factors in the SL(2) sector, meaning form factors 
with a composite operator from the above SL(2) sector. Throughout this paper, we will 
mostly focus on minimal form factors, i.e. on the case where we have as many external 
fields n as fields in the operator L.? Moreover, we will restrict ourselves to the planar 
limit; the employed methods are also applicable for finite N though. 

At the technical level, the unlimited number of covariant derivatives in the SL(2) 
sector implies that the form factor integrand at fixed loop order involves arbitrarily high 
powers of loop momenta in the numerator. This is a completely new feature as compared 
to amplitude computations or to the previous studies of form factors with operators in 
compact sectors. As demonstrated in this paper, on-shell methods can still be efficiently 
applied in this case. 


‘Note that the term magnon usually refers to a spin-wave. We use this notation due to the strong 
connection to spin chains and integrability, see [48] for a review. 

2We only consider L > 3 here. For L = 2, the remainder of the minimal form factor is a pure number 
as all dependence on the kinematics factors out. 


At one-loop order, we derive the full integrated form factor for all composite operators 
in the SL(2) sector. At two loops, we obtain the full integrand for an arbitrary magnon 
number using the unitarity method [49-51]. Solving the resulting integrals via traditional 
methods becomes more and more challenging for increasing magnon number. In this paper, 
we present the complete integrated result up to total spin three. We leave the solution for 
an arbitrary magnon number for future work. Besides universal IR divergences, the form 
factors contain UV divergences which allow us to extract the dilatation operator of the 
theory; our results (up to three magnons) match perfectly with the previously obtained 
data for the dilatation operator given in [52-54]. Despite the relatively small magnon 
numbers, we find very interesting patterns for the two-loop remainder function as well as 
for the finite terms at one-loop order. These observations indicate that form factors in 
N = 4 SYM theory have ‘hidden’ uniform maximal transcendentality for any operator, as 
will now be explained. 

The first pattern to be noted is that the remainder function contains a universal piece 
of maximal transcendentality four, which is identical to the BPS result [17]. This was 
first observed in the compact SU(2) sector and conjectured to be a generic feature of the 
theory in [21]. So far, further support for this conjecture was already found in the SU(2|3) 
sector studied in [28]. In the SL(2) sector, also lower transcendental functions occur for 
remainders of non-protected operators, i.e. combinations of polylogarithms and constants 
such as (3 and 7”, whose total degree does not add up to four. While this is expected from 
the SU(2) and SU(2|3) sectors [21, 28], a curious observation of the present paper is that 
the functions of lower transcendentality come with coefficients composed of (generalised) 


harmonic numbers or factors of ——, where m;,n; are magnon numbers as in (1.3). In 


particular, if also the latter are assigned an appropriate degree, the two-loop remainder 
is of uniform transcendentality four; similar structures can be seen in the one-loop form 
factor results. We refer to this assignment as hidden transcendentality.? 

In order to illustrate the above idea, let us give the following list of exemplary expres- 
sions with a (formal) degree of transcendentality k: 


Here, u; and v; represent kinematical variables, my, n; are magnon numbers as in (1.3), 
and the dimensional regularisation parameter £ is assigned a transcendentality degree —1, 
as normally done in the literature. The property of hidden uniform transcendentality can 
already be observed for the explicit expression for the two-loop dilatation operator given 
n [53]. The main observation made here is that this property appears to hold for the 
full form factor including functions with kinematical dependence, such as the remainder 
function. Note that in the case of the latter, we allow for generalisations of harmonic 
numbers or inverse powers of magnon numbers by inclusion of kinematical variables, such 
~ 3Note that the assignment of transcendentality to harmonic numbers is standard in the context of 


anomalous dimensions, see e.g. [55]. An assignment of transcendentality to similar fractions occurred in 
different contexts in the literature, see e.g. [56-58]. 


as . —. in the list above. A more detailed explanation as vvell as explicit examples 
will be given in section 4.3. 

The above leads us to the following conjecture: The full form factor for any given 
operator in N =4 SYM theory has hidden uniform transcendentality. 

Note that, in sectors beyond SL(2), the notion of magnon number might need to be 
suitably generalised, e.g. to the total number of excitations in the spin-chain picture for 
composite operators. Further note that the property conjectured above is not visible from 
the study of operators in compact subsectors as performed in [20, 21, 28], unless they are 
embedded into a larger non-compact sector (i.e. through inclusion of an arbitrary number 
of covariant derivatives). 

Further backing for the above conjecture comes from the relation between SL(2) op- 
erators and lightlike Wilson loops. It is known that twist-two operators (i.e. length-two 
operators in the SL(2) sector) in the large spin limit are related to lightlike cusped Wilson 
loops [59, 60], which have uniform transcendentality [61]. This is related to the fact that 
harmonic numbers S) become transcendental numbers of degree | for m — oo. Similarly, 
the generalisations of harmonic numbers become (poly)logarithms; e.g. .. -—— > 
— log (1 — it) for m — oo. Generic SL(2) operators are related to multi-edge lightlike 
Wilson loops. Since lightlike Wilson loops, which are dual to amplitudes, have uniform 
transcendentality in W = 4 SYM theory [62-65], this provides more support for our con- 
jecture. 

While the theory underlying the considerations of this paper is W = 4 SYM theory, 
we expect that the above insights may also teach us interesting lessons about similar 
calculations in QCD. In particular, one may wonder whether the maximal transcendentality 
principle [55, 66] extends to the notion of hidden transcendentality indicated above, see 
also [10]. We will further discuss this point in the outlook in section 5. 

The paper is organised as follows. In section 2, we give the tree-level form factors 
needed for the unitarity computations. Then, the full one-loop minimal form factor is 
obtained in section 3 and the one-loop dilatation operator and finite terms are extracted. In 
section 4, we construct the complete two-loop integrand and obtain the two-loop dilatation 
operator and the remainder functions for operators up to total magnon number three. We 
moreover discuss the interesting pattern and transcendentality property for the remainder 
functions. A conclusion and outlook is given in section 5. In appendix A, we give technical 
details on evaluating the one-loop form factor integrals. The complete remainder functions 
up to total magnon number three are given in terms of their densities in appendix B. 


2 Tree level 


In this section, we provide the tree-level form factors in the SL(2) sector which occur in 
the unitarity cuts in the subsequent sections. 


Minimal form factors. Within the SL(2) sector, a minimal form factor for a single- 
trace operator of length L contains L external on-shell scalar states. While the state 
can be conveniently expressed in terms of oscillators acting on a (spin chain) vacuum 


(671, the minimal form factor is simply obtained by replacing these oscillators by spinor 


helicity variables (A, A) and the SU(4) Grafmann spinors 74 [19], see also 168, 69]. For 
3 


the SL(2) sector, each scalar X = 34 at position i contributes a factor nnf and each 
covariant derivative D* contributes a factor pi = pf ə For explicitness, but without loss 
of generality, we choose T, = (ay), such that Dt = DH and pt = p! = MA!, when 
required. The minimal form factor can thus be written as 

L L ( + eo 


L 
F o (En) X pe 
1=1 


: The j Mkt j > (2.1) 
kel \j=1 ` 


where the sum over k is due to the cyclicity of the colour-ordered form factor, which in 
the minimal case is essentially given by the manifestly cyclically invariant vertex of the 
operator. Having fixed the type of scalar X, a given operator or minimal form factor, 
respectively, is hence characterised by the derivative configuration 7 = (N1,...,NL). 


Next-to-minimal form factors. As important building blocks for the two-loop com- 
putation via unitarity, we will also need the next-to-minimal tree-level form factors, which 
have L + 1 external states. These can be easily computed via Feynman diagrams. The re- 
sults coincide with the corresponding tree-level expressions for general operators obtained 
via the twistor action in [26]. 

At MHV level, there are two different configurations of external states to be considered 
for the next-to-minimal form factors. The first configuration consists of L—1 scalar external 
states and two fermionic external states. The corresponding next-to-minimal form factors 
can be obtained from the minimal form factors by replacing X; = pa > we a — uduş 41 
in the external states of the minimal form factor (2.1) for each į = 1,..., L according to 


göyə a y —— y, PD" hy 
Volt 


(ii+1) ml (nj — m)l 


(neni — MM) > (2.2) 
m=0 
where we also have to relabel pr > Pia for l = i+ 1,..., L and include the momentum 
PL+1 in the momentum-conserving delta function. 

The second configuration consists of L scalar external states and one positive-helicity 
gluon external states. For this configuration, two qualitatively different cases occur, which 
are distinguished via the SU(2) indices, i.e. the indices & of A. In the first case, all SU(2) 
indices remain on the scalars. This process can only be attributed to a pair of neighbouring 
scalars, not to a single scalar, as gauge invariance is achieved by summing two Feynman 
diagrams: one with the left scalar emitting the gluon to the right and the other with the 
right scalar emitting the gluon to the left such that the gluon appears between the two 
scalars in the colour-ordered expression. We can obtain the contribution of this case from 
the minimal form factors (2.1) by replacing X;Xj41 > X; gözü in the external states 


for i = 1,..., L according to 
: + : = pk E 
məx 3,4 (Pin B a ə — 52 (pz yə ney” urina, (2.3) 


“Note that we are suppressing a factor of the imaginary unit i from the (covariant) derivative. 


where we also have to relabel pr > m”: for l =1+2,...,L and include the momentum 
pL+1 in the momentum-conserving delta function. In the second case, at least one SU(2) 
index from scalar 7 lands on the gluon. Two Feynman diagrams contribute to this case: the 
emission of a gluon by scalar 2, and the insertion of a polarisation vector into the gauge- 
field part of the covariant derivative acting on this scalar. Their sum is gauge invariant, 
and, as they are only involving the scalar ?, we can write this contribution by replacing 


Xi > XG + g7 Xizi in (2.1) for i = 1,..., L according to 


zi nj ə = 
tiyə un Su MAH : 1 gyre (pu) i 3,4 
nj! (ii + 1) £ m(nj —m)! (m— 1)! 7: 
=i . (2.4) 
n ge As — 
AAi L1 (pz )™ i (Pia) r nf 
E s “07, 431 £., ...80 ED Ge os 
(1-1) <= m(m— 1)! (nj -m)l -— 
where we again relabel pE > Pia for l =i+1,..., L and include the momentum pz+1 in 


the momentum-conserving delta function. 

We acquire the next-to-minimal MHV form factor by applying the replacements rules 
(2.2)—(2.4) to all fields in the operator and summing the resulting expressions. The next- 
to-minimal NMHV form factor can be obtained from its MHV counterpart via conjugation, 
as described in [20] in our conventions. 


3 One loop 


In this section, we calculate the one-loop correction to the minimal form factors in the 
SL(2) sector and extract from them the one-loop dilatation operator and the finite part. 


3.1 Unitarity construction 


Let us write the loop expansion of the minimal form factor in the following form: 


Fos (1 +5 Pero) : (3.1) 
=i 
where : 
GeuNe B 
2 ar > (Are By (3.2) 


is the effective planar coupling constant. For operators such as BPS operators or the 
Konishi primary, which are eigenstates under renormalisation, Z is simply the ratio of 
the £-loop and tree-level form factor. However, this is not true in general due to operator 
mixing. In particular, the loop corrections to vanishing tree-level form factors can be 
non-vanishing. It is therefore important to promote Z to an operator that acts on the 
tree-level form factor FY ) and can create different tree-level form factors from it. Examples 
in the SU(2) sector and for general operators can be found in [21] and [19], respectively. 
In the planar limit, connected £-loop interactions can maximally involve £ + 1 neigh- 
bouring fields in the colour-ordered form factor at a time. Hence, Z can be written as 


Pl 


P2 


Figure 1: The double cut of one-loop form factor in the s¡>-channel. 


a sum over connected interaction densities each of which involves £ + 1 adjacent external 
legs plus disconnected products of lower-loop interaction densities. At one-loop order, the 


maximal interaction range is two, and we can write" 


L 
25377 (3.3) 
4=1 


(1) 


where 1;;\{ acts on the external fields ? and 7 + 1 and cyclic identification ¿+ L ~ i is 


(1) 


understood. Ij}; can be given in an operatorial form as 


co +ym + ma ni na 
oe m y (E yr ) +1 ( ) ( (3.4) 


! ma! + + ” 
ni,mi=0 ma 2 Op; pi =0 Op; +1 pi, =0 


where (J;)njino” is the matrix element of the interaction density that corresponds to chang- 
ing the derivative configuration at momenta p;, p;+1 in the tree-level form factor (2.1) from 
ni, na to M1, mə, respectively.” Here, (È) should be understood as acting to the right 


such that yn Ər F(p) 
p 


p=0 


(3.5) 


ı 
p=0 
where the action is understood to be on the polynomial in P}, not on the momentum- 


conserving delta function. 


1 : ! 
))m1;m2 can be computed using the uni- 


One-loop interaction. The loop correction cr 
tarity method. We will be comparatively brief here; see e.g. (19—21, 37] for details of similar 
calculations and 121) for details on our conventions. As only interactions between two adja- 
cent fields occur at one-loop order in the planar theory, it is enough to consider the double 


cut shown in Fig. 1. The cut integrand takes the form 
farsa, l2) dim, Am, FS oe ls, p3, «PL; E 27:75 (3.6) 


[farsante ar dali y ( 2 a FS (pı, -pu:4). 


nı! məl (11) (22) Öp? p1=0 Ops p2=0 


5We use the italic J for the density and the calligraphic Z for its sum. 
SNote that for compactness we use only the index i instead of i,i + 1 for the symbol I if the range of 
the interaction is obvious, e.g. due to explicit specification of the spin numbers. 


where the minimal tree-level form factor is given in (2.1) and the four-point amplitude is 
given by the standard MHV expression. The two-particle phase-space measure is denoted 
by dLIPS. 

The full integrand can be obtained from the cut integrand by adding the cut propa- 
gators and taking the cut momenta İrə to be off-shell. Explicitly, the square bracket part 
on the right hand side of (3.6), now written for general i, lifts to the one-loop integral as 


Pi 
TOY s. > ("Gy 
(Li hmi nə = “Suzi ni! nəl 
Piri (3.7) 


¡(0% p2)? : TH (Deye 
i mln! J ir? PBG — pi)?” 


where lo = pi + pi41 — li. The 't Hooft mass is denoted by u, the spacetime dimension is 


= =$ 


D = 4 — 2e and yg is the Euler-Mascheroni constant. This expression agrees with the one 
already found in [19]. 

For given small magnon numbers mi and mə, it is straightforward to perform standard 
Passarino-Veltmann (PV) reduction [70] and integration-by-parts (IBP) reduction [71, 72], 
which reduces the integral to the scalar bubble integral (A.10). The integral (3.7) evaluates 
to a polynomial in terms of the lightcone components of the two external momenta as 


ni+na +ym + yma 
1 1 (pm (pu) 
m 2. Oe mr (3.8) 
m1,m2=0 1. 2" 


my+m2=n1+n2 
As shown in appendix A, the coefficients in the above expression can even be evaluated 
to a closed form in the generic case: 


Pi 
¿Oyumoz mamar (1-2) = 4 
200 HOLD +g = 1 ES 1) m 


Ru AE +m +n 31) 


J! (mi — j)! (ma — j)! 


min(m2,n1)—j R z 
y m —? mə —) 
nı— j-r r 


r=max(0,m2—n2) 
T(2 + ng — mə + 2r — 2)7(? + ny + mg — 2j — 2r — 1) 
(3.9) 
where the diagram denotes the standard scalar bubble integral given in (A.10). We have 


j=0 


checked in numerous special cases that the general result agrees with the one obtained 
using PV and IBP reduction. 

The result (3.9) contains the full e-expansion, which, for example, also includes the ra- 
tional contribution.’ This will be important for computing the two-loop dilatation operator 
and remainder function. 


"The cut-constructible part has been computed in [19]. 


3.2 Dilatation operator 


The anomalous dilatation operator is defined via the renormalisation factor Z according 
to 


ôD = == log Z = reg? A log Z. (3.10) 
du Og 
Expanding in powers of the coupling constant, we have 
22305 z— Y gzo, ZO =1. (3.11) 
(=1 (=0 


At one-loop order, the renormalised interaction (indicated by the underscore), 
TU = 74 4 zü. (3.12) 


is UV finite by definition of the renormalisation matrix Z(), As TÚ), we can expand ZA) 
in terms of densities as 


L 
zülm T (3.13) 


(1) 


where we expand Z;;,, in terms of its matrix elements in complete analogy to (3.4): 


oo +\m + \M2 nı n2 
A= 5 JE POT Wr)” ( ) ( ) | (3.14) 


! ! T + 
ni,mi=0 1 2 Pi / př =0 NOPi41 Pi =0 


By definition, 20 has to cancel the UV divergence of J;;,1. In terms of matrix elements, 


m 


where we have written the UV divergence as the total divergence minus the IR divergence. 


this reads 


(20) yma ma = (1P yma me 


a n1,n2 n1,n2 


uv Lo 2. 


n1,n2 n1,n2 


(3.15) 


> 
+-pole 


The former is then given as the remaining +-pole, which we extract as indicated by the 
vertical bar. The IR divergent part is universal and given by (see e.g. 1191) 


Pi —_E 
cr! yayma 5 = “0003 —T — m” (- Fm ) 0.0 + Ole ), 
Pi+1 


(3.16) 
where the scalar triangle integral is given in (A.12). 
From eguation (3.10), we can read off the one-loop dilatation operator as 
1 1 1 m m 1 m m 
.—. 2. = wA a (a! Dyna = 2e(Z/ rie (3.17) 


where we have expanded the one-loop dilatation operator in terms of densities and matrix 
elements in complete analogy to (3.13) and (3.14), respectively. We have checked in explicit 
cases that the result from the above expression agrees with the well-known expression for 
the one-loop dilatation operator [67]: 


n1,n2 2 — mə ni+n2 


(DO ymasma =? Ea gm2 — zun qn sof mEn ) . (3.18) 
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Here, we employ the harmonic numbers defined by 


n 
1 
002 r (3.19) 
k=1 
as well as the Heaviside function 
1 forn>m 
Onm = , 3.20 
7 | 0 forn <m; ( ) 


see also [53] for the notation used here. 


3.3 Finite terms 


We would like to better understand the structure of the complete one-loop form factor 
obtained above. After extracting the dilatation operator as the 1/e-pole of the interaction 
symbol Z, let us now study the finite part of the form factor in the e-expansion. We define 
the non-trivial finite part as the O(<%) contributions in the e-expansion of (3.9) after the 
UV and IR divergences are subtracted using (3.15) and (3.16): 


Fin (Cae =) = (1P yemə — (¡Dynamo 4 (gyms, (3.21) 


n1,n2 n1,n2 Ting 


This yields in the special case mə = 0° 


mim magma |_1 2 3 İİ 
.—.:125:5— [+ 


1 m m 
+ a” log ( — o (8.22) 
The general case can be obtained from this special case via the folloving symmetry relation 
for Hö: 
my (Tp? msde + ma (Imama = (m + YOM ETA, + (m2 +) a (3.23) 


which is a direct consequence of momentum conservation and leads to an identical equation 
for the finite part. As a result, we find 


Fin Gone = (3.24) 


NIMA 


ng n2—) 


mı! ma! 1)) ng! (nı +3) Minds ent 

my! ng! 2, D m (na — j — k)! (mı — k)! (ma — nə +1)! -. (U; —. ) i 
where l = j + k. 

We can see that the finite terms in (3.22) and (3.24) have a simple structure similar 
to the one of the dilatation operator in (3.18). In particular, they have uniform max- 
imal transcendentality two if the harmonic numbers and inverse powers of differences 
of magnon (177 are 2 a degree of 2... according to die .. 
su 


eee 


8Note that we could have absorbed the logarithm term by adding (2? mm ( — “əmi instead of 


n1.n2 
1 
(em, 
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4 Two loops 


In this section, we continue the study of SL(2) form factors at two loops. We will see 
that the unitarity method allows us to obtain the full integrand. We solve for the explicit 
dilatation operator and the remainder functions up to spin three. The two-loop dilatation 
operator and remainders show interesting patterns and suggest a hidden uniform maximal 
transcendentality property. 


4.1 Unitarity construction 


We first compute the two-loop integrand via the unitarity method. At two-loop order, 
connected interactions involve at most three fields of the composite operator, which are 
adjacent at the planar level. The two-loop corrections can be expressed as 


L 
mas 57 5: +5 È Ta 1 (4.1) 


i=1 


(2) 


where f;;, , and Z, ria ¡+2 denote connected interactions of range two and three, respectively. 
We expand the interaction densities in terms of matrix elements as 


+ 
rə) = Y a A ( 2 ). (4.2) 
11+1 7 i /nı,mə my! məl Op; py =0 Opes he 
10) ee Y (IP y mm me (py y” ! (pia) (phəy"" 
ii+li+2 a 2 7N1,N2,N3 my! məl məl 
mn) 
Đ nı ð n2 ö n3 
rl, Ca Ce ), SE 
Pi / py zo \9Pi41/ pi, =0 Pia P;4270 


The one-loop correction (I CO) ? has been given in the previous section. We apply 


the same unitarity technique to compute mö and tm ia the corresponding 
non-trivial double and triple cuts are shown in Fig. 2. In the following, we briefly outline 
the major steps and present the full integrand. More technical details on this kind of 


calculation may be found in [21]. 


Double cuts. The simplest cuts are the two-particle cuts in the two-particle channels, 
with the one-loop amplitude or form factor on one side of the cut. They are shown in 
Fig. 2a and 2b, respectively. From them, we can fix the integrand contributions from the 
following topologies as 


ly [t+ nı I+ na It n3 
2 
in = +1 aa — 2 ) ( 3 ) 
l3 s 


i+2 
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(a) First double cut in the s12-channel. (b) Second double cut in the s12-channel. 


q P1 


DA P3 


(c) Triple cut in the s123-channel. (d) Triple cut in the sış-channel. 


Figure 2: Unitary cuts of the two-loop form factor .. 
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where n1, mə, ng are the numbers of derivatives occurring in the corresponding minimal tree- 
level form factor and thus in the integrand. Moreover, we have abbreviated lg = p¿+Pi+1—l1 


2) 2) (2) 
o 


in IO) and la = pi +Pi+1 +Pi+2— lı — l3 in I as well as I; . Note that these contributions 


can be equally obtained via consecutive double cuts. 


Triple cut in the three-particle channel. Next, we perform the triple cut in the 
three-particle channel as shown in Fig. 2c. Both (4.5) and (4.6) contribute to this cut. The 
new contribution determined by this cut is 
lı I+ nı i n2 It n3 
2 
ñas => dd 27 1 ) (2 05) 
l3 ga nı: 


(4.7) 


Triple cut in the two-particle channel. The most complicated cut is the triple cut 
in the two-particle channel as shown in Fig. 2d. All previous integrands contribute to this 
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cut. The new integrand contributions are 


A i+1 h i+1 
(P) — ( + +ph j= + 
¿¿EJn1n2n3 — Dizi T Pirə Pj415le ls 
l3 i+2 la i+2 
i+1 


lı nı 1 (pf” yı Xx05000 
= : : + + 4 WPi) (4) Uz 3 7 
> = Si+ti+2lli + ) 3 m(ni=m)!(m=1)! nəl nal (4.8) 
m y EDO DA DIN Gy" Xə 1 öy" E DA (yn ED 
=L ee Bs. S a m ni! (na—m)!(m—1)! n3! 
p$ Lehn gön aye gön 
mər nı! nəl (m — 1)! (nz — m)! 
and 
(2) 7 i 2 i 
(Ti F)nuynana = -< (py + py) = —T piya l2 
l3 itl l3 gi 
l i nı Jimm Amei oe an 
1 (yn (A (A ia)” 
= aq (İ +I 1. 3 i T 
l3 . +1 (1) D) SAA na! nəl (4.9) 


as vvell as 
(19) x aa tf: a O 
4,G mi mə in "= ii+l i+1 3 E m (nı — m)! m, — 1)! nəl 
boty yy SD" ae UP? qiy" qiy zan 
“m. Lom m! (m-1) (nz — m)l mil nəl l 
Here, we have lg = pjy1 + pi+2 — li — l3 in i and l2 = pi Epi — li — lg in rg as well as 
ır. 


We see that the structures of the loop momenta in the integrand are directly inherited 
from the corresponding structures in the tree-level form factors (2.1)-(2.3). 


Result. Given the above building blocks, the interaction densities can be defined as 


12), = 141, m 
COE +12 +14 41, | 


2, 2, 
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They can be expanded in terms of matrix elements as 


mə 


Umm = Y MO meo LO (Pi)? 


$ 11542 my! mo! 
(4.12) 


CP maura = E (Pan E (elə Gah 
— E - mil məl məl 
where mi + ma = mi + na and mai + ma + m3 = mi + na + na, respectively. 

Let us explain the special combination of building blocks in (4.11). We should first note 
that there is an ambiguity in grouping and distributing the terms with integrals involving 
only two external momenta. In (4.11), a choice is made such that the terms in T . and 1 i 
are interpreted as contributions to the interaction density of range two acting on the sites 
i and ¿+1, while the remaining terms are understood as interactions of range three acting 
on the sites 2, ¿+1 andi+2. While this assignment is irrelevant when summing over all 
densities to obtain Z2), it ensures the finiteness of the remainder density obtained from 
the corresponding interaction densities. We will come back to this point in subsection 4.3, 
see the discussion after (4.25). 

As in the one-loop case, the integrands (4.4)—(4.10) contain uncontracted loop mo- 
menta. Using PV reduction, they can be expressed in terms of Lorentz scalar integrands, 
which contain loop momenta in the numerator only in the form of irreducible scalar prod- 
ucts. The resulting integrals can be reduced to master integrals via IBP identities, using 
e.g. LiteRed [73]. The required master integrals are all known analytically and can be 
found in [74]. This procedure can be efficiently applied for small numbers of magnons; 
we have computed explicit results for all cases up to a total magnon number three, i.e. 
¿ni x 3. While in principle applicable to higher numbers of magnons, the computation 
time for the employed method significantly increases from four magnons on. 


4.2 Dilatation operator 


The two-loop renormalised form factor is given by 


10 — 70) 47020 4 Ze), (4.13) 


where the two-loop renormalisation constant has to cancel the UV divergence and is given 
as 


ial 


295 (Zaa +5 = 5 29,29). (4.14) 

i=1 = =i+2 
The first term inside the summation above is due to the connected interactions, while the 
last term accounts for all disconnected products of one-loop interactions in (4.1). Unlike 
in the one-loop case where the IR and UV divergences can be easily separated, at two-loop 
order the IR and UV divergences are entangled with each other. We can subtract the 
(universal) IR divergences via the BDS ansatz [75, 76], as explained in detail in the next 


subsection. The two-loop renormalisation constant density Ze 42 can be obtained by 


imposing the BDS remainder density RO) defined in (4.25) to be finite. 


ii+li+2 
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In terms of the renormalisation matrix Z, the form of the two-loop dilatation operator 
follows from the definition (3.10): 


1 
DO = de (zƏ = (Ey) . (4.15) 
It can be expanded in terms of densities as 
= 
0 (4.16) 
i=1 


As usual, the densities are expressed via matrix elements in complete analogy to (4.3). Via 
the above relations, we may extract the matrix elements of the two-loop dilatation operator 
in analogy to the one-loop case. Note that the square in (4.15) has to be evaluated as a 
square of matrices. As an example, consider the following matrix element of the dilatation 
operator density, which is given in terms of the (non-vanishing) matrix elements of Z() 
and ZÜ) as: 


2 


al(ss m) “909 | = a 


where the boldface number indicates the magnon number at the single internal leg. The 


1 
OPES EE EEN 


matrix product will also be used later in the definition of the remainder function, see (4.25) 
for a pictorial representation. Instead of giving explicit results for all matrix elements up 
to three magnons, we may compare the obtained expressions to the previous literature, 
where the two-loop dilatation operator in the SL(2) sector can be found in closed form. 


Comparison with Feynman calculation. Notably, the full two-loop dilatation opera- 
tor in the SL(2) sector of M = 4 SYM theory has been obtained by Belitsky, Korchemsky 
and Mueller in [53]. The named authors provide an explicit expression in closed and rela- 
tively compact form. Our computation of the two-loop dilatation operator for states of up 
to three magnons shows perfect agreement with that operator if we evaluate both expres- 
sions on composite operators. For convenience, we display the respective expression in the 
following.’ 

In contrast to our expansion (4.16), the authors of [53] write the two-loop dilatation 
operator in terms of a range-two density and a range-three density:!° 


L 
PO = (Dia + 2) - (4.18) 
i=1 


“Note that the expression printed here differs from the expression given in [53] by a conventional overall 
factor of —2 and by some conventional factorials of magnon numbers; we make these factorials explicit in 
the definitions of the matrix elements, which are defined in complete analogy to (4.2) and (4.3). We also 
insert the correct factor of 2 in front of so se in (4.19) compared to the expression given in the appendix 
of (53). We thank Grisha Korchemsky for communication and sharing related computer algebra files with 
us. 


10Note that one may include the range-two density DO into the definition of the range-three density 


(2) 


Oe ced by formally adding identity legs to 200 
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The range-two density reads 


~ Ün m: m m 
(OP = 2/2509 si) + s) oman + — ama smn [SO +380) 
+ (58) — s) (s) +350) — asl) mom) | HAA, (639) 
and the range-three density reads (for mə = 0) 
u ¿mi+ma+ma y 
(2) ymi,mə ma — _ n1+0+n3 n3zm3 1 0) — qü). qü) 
(9; Inia 2 (nı = mı) (ns = mg) q (= -mı + ə 7 ro ml 
1 nyen 
+ Onimi Ca mi + Se m3 sm) a { mem > (4.20) 


as well as (for mə £ 0) 


n2 n2-j1 n2-jn2-j=k1 


( 90 ym ma ms 27 1)/nəl 
-— Ima! 2 2 Ja 2 = ! k3! (na — £)! 
711: 71217) 7) : 
1: 102. 3! j1=0 j3=0 kı=0 k3 ii! 3 1:3: ( 2 ) 


(ni + ji)! (n3 + 53)! = (2) m1 —k1,m3—n9+8,m3=k3 
(mi — k1)! (ma — nə + £)! (m3 — k3)! (D; trios tia » (4.21) 


where j = jı + j3 and £ = j + kı + k3. Note that the match with the two-loop dilatation 
operator of [53] is only at the level of the composite operators, i.e. after we sum over 
all matrix elements of the dilatation operator density that occur for a given (cyclically 
symmetric) single-trace operator. The matrix elements of the dilatation operator density 
do not directly match with the expressions of [53], which is of course also not required. To 
make this clear, we have marked the matrix elements of [53] with tildes. 

We notice that the matrix elements (4.19) and (4.20) of the two-loop dilatation operator 
are of uniform transcendentality three if we formally assign the transcendentality k to the 
- for j = 1, 2,3. 
With the same assignment, the matrix elements (3.18) of the one-loop dilatation operator 


harmonic number SÜ”) and the transcendentality one to the fraction = 
2 


are formally of uniform transcendentality one. Let us refer to this property as hidden 
uniform transcendentality. Note that the assignment of transcendentality to harmonic 
numbers is standard in the context of anomalous dimensions, see e.g. [55]. An assignment 
of transcendentality to similar fractions also occurred in other contexts in the literature 
[56-58]. 


Comparison with algebraic construction. In [54], Zwiebel introduced an algebraic 
construction that allows to recursively compute the two-loop dilatation operator within the 
SL(2) sector. We have implemented this loop-recursion and we compared the spectrum of 
the resulting operator to the above findings for the M = 1,2,3 excitation sectors. In all 
cases, we found perfect agreement. The operatorial form of the two-loop dilatation operator 
obtained by the above recursion is rather lengthy and we thus refrain from displaying it 
here. We note, however, that the recursion in [54] actually applies to the larger PSU(1, 112) 
sector of W = 4 SYM theory, which makes it interesting for future extensions of our work. 
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Another advantage of this algebraic method is that it allows to construct a dilatation 


operator that is manifestly hermitian.!! 


4.3 Remainders 


The (two-loop) BDS remainder [75, 76] in the refined case of renormalised form factors [21] 
takes the form 


RO = Te) - > (LV) - f(e)Z (gə) + Of), (4.22) 


where 
FP) = -2 — 2e = 2422. (4.23) 


Similar to the example discussed in (4.17), the product of the one-loop density with itself 
has to be understood as a matrix product. For the two-loop remainder 


L 
2 
RO) = Hn (4.24) 


this is captured by the folloving pictorial expression 


2 
Rİ 1 i+2 — 


K NG 77. 
+-| = - 1 (560) 1-5 65) 
2 CON : ? ce 


2 
The remainder density Re i la isə 18 expressed in terms of its matrix elements (RÍ aramı in 


complete analogy to (4.3). We determine 2: ¡+2 from the requirement that the remainder 
density RG), ¿+9 is finite. This is possible only if momentum-dependent 1-poles are absent 


in Rİ . i+2 
(4.11) of range and to average the range-two densities as in (4.25). 

We have obtained explicit expressions of the remainder function (4.25) for all configura- 
tions up to magnon number three, i.e. (RO yma, 1 with nı +n9+N3 = Mi] +ma+mg < 3. 


The densities can be written as linear combinations of 14 functions, which are partially 
related by relabelling the kinematical variables. The latter also enter in the form of 


before adding z% a 42, for which it is crucial to have the correct assignment 


Sii+l Sj4+1i+2 Sii+2 
YU = VE) 0, — ——, wy ——, (4.26) 
Sii+li+2 Sii+li+2 Sii+li+2 
where 
Sii+li+2 = Sii41 + Sit dita + Sii42- (4.27) 


“Note that in perturbation theory loop corrections to the dilatation operator do not have to be hermitian, 
see e.g. [77, 78]. 
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Concretely, we can write 


255 :51457 
kol 
where İ(£, ay, na “İk are rational functions of ui, vi, w; and 
2 2 
7 2 2 2 T i+] i TU iyi 
Y, = G ES RN as s AE qe o 2), — los ( AH), 
= 
log ($) log (— =) + — log (22) log ( — m... s) 
1 z? 
Li. (1 — v,) — a log? (u;) + log(us) log(v;) — TE 
. 1 2 72 
Lig(1 — uj) — j log” (vi) + log(u;) log(v;) — m 


2 
- © ,log(ui), log(vi), log ( = -. 1) ; 
(4.29) 
In the vector V;, the components are ordered according to their (decreasing) degree of 
transcendentality. The first component of V; has transcendentality four and is given by!” 


1) — log (45 we 1) — Lig(1 — wi) 


1 
a Log (vz) —— 7 og” (1— u) 


= (RO) XXX = —Li,(1 — si) — Lig(u;) + Li, (5 


— log (ui) |Lig( 42) + Lig( — 22) + Li: (51) — 


“Vi( 42) Di (r) + Lia (u) flog (53) log (vi) + log? zi) 
+ log (us) — E log? (us) log? (vi) — $ log? (1 — ni) log (us) log (£) 


1 1 
-zlog(1— ui) log? (uj) log (vi) — ; log” (uj) log (w;) 


— m sl a 
= m [og (us) log (13) + 5 log? (13%) — 5 log? (us) 


Se toate u 200 A 001 (4.30) 


It is the density of the BPS remainder first obtained in [17]. Moreover, 


w2 


r 5 log (1 — ni) log (vi) log (q 2i) 


= iat — at) — log (ui) Lio (3 


| 
ost a. 7 los (ui) log (vi) log (wi) — = los? (wi) + (ug + vi) 


Wi 


1 
— Lig (1 — vi) + Lig (ui) — 5 log” (vi) log (E + =r’ log (Y E) 


1 2 Si i4+1i42 
m 6" log ( — —. 


(4.31) 


12 Apart from classical polylogarithms, this expression contains the Goncharov polylogarithm 
G ({1 — Ui, 1— Ui, 1,0} , 02). 
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is the (up to relabelling) unique transcendentality-three contribution that occurs for re- 
mainders in the SU(2) sector built from the complex scalars X and Y [21]. Here and 
below, |aeg-k denotes the projection to the part of transcendentality degree k. Note that 
at transcendentality two, besides the (up to relabelling) two functions encountered in the 
SU(2) sector 


1,5 > > > = 
(ROEE = gl Vile + [Vaz + [i]s + [Velo — Wiho, 
: La La B (4.32) 
(Röya ə” 3 IV Me, 


only (up to relabelling) one further transcendentality-two function is required to express 
the transcendentality-two part of the remainder in the SL(2) sector. A similar relation 
holds for the functions at transcendentality one. Hence, compared to the SU(2) sector, an 
extended basis of functions with coefficients that contain ratios of Mandelstam variables 
is required. This goes also beyond the SU(2|3) case, where in addition to the remainder 
densities in the SU(2) sector with rational coefficients only the constants vr? and Cs were 
encountered in [28]. The additional functions are still rather simple though. 
One trivial symmetry of the remainder density is what is commonly called parity: 
RO mamma — (RO) əzənə ) 4.33 
( EU n2,n3 = ( 2 7-5 Si itO Siti it Uii ( ) 
which is a simple relabelling.!? This symmetry is already visible at the level of the integrand 
of the two-loop interaction density. Moreover, as a consequence of momentum conservation, 
the two-loop remainder obeys a relation analogous to (3.23) at one loop. It reads 
m (RP gema RO a” ma Ry Yana (4.34) 


= (nı + Dry, mams (ii Dre mamama + (ng + 1)(RO jr mama 


i /nyti,ng,ng i /nin2+l,n3 i J/ni,n2,n3+1 * 


Using this identity, we can reconstruct (RÜ yy məmə from the knowledge of (RÜ) mama 


in complete analogy to (4.21): 


m mol mal <2 n2—j1 na—j na—j—k1 
(RO ym: məmə — 1 2 PER EQ hs 3 3 3 2 
mmama ny! gl ng! 


777 


Ling! 
ji! ate (n2 — £)! 


HET res a ee, (4.35) 


where j = jı + j3 and £ = + kı + ks. Thus, it is sufficient to look at (RÜ) jamama with 


2 ny 0,n3 
m1,m2,mM3 


nı > ng, and mı > mə if nı = ng. The corresponding coefficients (€?) 0 ns 


are given in 
appendix B. 
Inspecting the explicit expressions in the appendix, we observe several interesting 
patterns. Firstly, the highest transcendental part is universal: 
Ghee =e On Ons (4.36) 


N1,N2,N3 n3 ? 


13Notice that the entries of V; organise into singlets and doublets under parity. 
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which is to say 
(RI 


2 n1,n2,n3 


eae |g = S! Oa One (Ri VERE (4.37) 
eg- 


This property was already observed in [21] for the SU(2) sector and conjectured to extend 
to all operators. Secondly, also the coefficients of the transcendentality-three functions [V;] 2 


> 


to [V;]5 follow some interesting patterns. We analyse the following expressions as examples 


below 
(2) ym,0,0 = m 1 1- wi \* | p(2)\xyx am T? _ SiHLi+2 
(R; 1500 deg-3 Xə w; ) (R; Vee dəə ə ( 1) 6 log ( m... 
= a 
(Vile [Vila 
(4.38) 
ara] = + (HE tog (- #2) (4.39) 
2 m,0,0 deg-3 = m 6 g m ? $ 
a 
Vila 
(2)\0,0,m — (-l)” up) xvx 4.40 
(Ri 7 0,0 deg-3 — m w” (Rì Wee deg-3 ( z ) 
Vile 
(2) m—n,n,0 = (2) m,0,0 (2) 0,n,0 4.41 
(Ri 7: deg-3 Ön (Eu A. de rs i 7: daşa ( s ) 


Similar patterns are encountered also in the coefficients of the transcendentality-two func- 


tions; e.g. 
(1) 2 (1) 2 
(2) \0,m,0 Sm Uj Sii+1 T Sm 1 T 
E ae = — | log (=) log ( — E) + — -=—=-—=|(-1)—. (4.42 
al, ee (o CH) do (— 25) + =) +( dl m JE. (442) 
5 
[Vile [Vi]1o 


We see that the coefficients of transcendental functions contain harmonic numbers and 
1 
mj=nj 


in the one- and two-loop dilatation operator (3.18), (4.19) and (4.20) and the one-loop 


factors +, where 4 = for some j. These are precisely the building blocks occurring 


finite terms (3.22). Moreover, we have generalisations of harmonic numbers which include 


ak 
ratios of Mandelstam variables: 3777 LL =) . For1— u; = w;, the latter reduces to SS), 


The given matrix elements of the remainder function density have a formal transcen- 
dentality degree four if we attribute a degree to the above coefficients according to the 
following rules: 


uy mi — Ni Mj — Nj Wi 


= — u,)” 1 1 i 
SO >l, 25 —I. —— 6 (4.43) 
a 


This generalises the property of hidden maximal transcendentality observed for the two- 
loop dilatation operator in the previous section to the case of functions of (ratios of) 
Mandelstam variables. In the limit of large spin, we have 


Nik 
xu... 


lim 
m— oo m— oo Wi 
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i.e. the result is of uniform maximal transcendentality in the usual sense. 

It would be very interesting to understand the full pattern of the remainders, in par- 
ticular to obtain a closed expression in analogy to (4.19)—(4.21) for the dilatation operator. 
In order to guess such an expression, we would need more data for higher numbers of 
magnons. We leave this point for future work. 


5 Conclusion and outlook 


In this paper, we have computed the minimal form factors of composite operators in the 
SL(2) sector of M = 4 SYM theory up to two-loop order. In particular, we have determined 
their finite remainders as well as the dilatation operator. The employed on-shell unitarity 
methods have once more proven to be an efficient tool in this context. Let us summarise 


the main results in some more detail: 


e At one-loop order, we have calculated the full minimal form factor to all orders in the 
dimensional regularisation parameter. Extracting the one-loop dilatation operator 
from this data, we found perfect agreement with the existing literature. Moreover, 
we have extracted the previously unknown finite terms. 


e At two loops, we have determined the two-loop integrand for states with an arbitrary 
number of magnons and the complete form factor remainder and dilatation operator 
for states with up to three magnons, i.e. with up to three derivative insertions into 
the composite operator. Also here, the obtained dilatation operator matches with 
the closed expression that was determined previously. 


e While the remainders contain functions of lower transcendentality, we have observed 
that their coefficients follow some interesting patterns. These support the conjecture 
of a hidden uniform maximal transcendentality property, cf. section 1. 


In order to confirm the above conjecture, the next logical step is to consider further 
examples of two-loop form factors. In particular, an important goal is to extend the above 
results to all composite operators in the SL(2) sector, i.e. to composite operators with 
more than three magnons. This requires a more efficient approach to the PV and IBP 
reduction of the general integrals (4.4)-(4.10). Although it would be desirable to solve 
the integrals for arbitrary magnon numbers, this is not strictly necessary. As the hidden 
uniform maximal transcendentality property conjectured above significantly constrains the 
dependence on the magnon numbers, we can fix the coefficients of an appropriate ansatz 
from explicit results with sufficiently high but finite magnon numbers. In order to better 
understand the transcendentality patterns, it may then be useful to consider the inverse 
Mellin transform of the remainder results.“ 

It would also be interesting to consider larger subsectors closed under renormalisation, 
such as the PSU(1,1|2) sector. The two-loop dilatation operator in this subsector can 
be computed via the iterative construction given in [54], which should provide a useful 


“For an example in the context of Wilson coefficients see e.g. [58]. 
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guideline. Finally, it would of course be very interesting to calculate the two-loop form 
factor remainder and the two-loop dilatation operator in the full theory, which is currently 
unknown. We believe that the study of SL(2) form factors, with their new feature of 
non-compactness, provides an essential step in this direction. 

While VV = 4 SYM theory is interesting in its own right, methods developed in the 
context of this model have often lead to broad applications, such as in QCD or pure Yang- 
Mills theory. Apart from this, results from M = 4 SYM theory are also useful building 
blocks in QCD. Notably, based on the so-called maximal transcendentality principle [55], 
the three-loop universal anomalous dimension was initially extracted [66] as the leading 
transcendentality contribution of the three-loop non-singlet splitting functions of QCD 
[79].1% Further evidence comes from the study of the three-point two-loop form factors of 
the stress-tensor supermultiplet in M = 4 SYM theory [10], where it was found that the 
two-loop remainder of W = 4 SYM theory matches exactly the leading transcendental part 
of the two-loop Higgs-to-three-gluon amplitude in QCD [82]. This generalises the maximal 
transcendentality principle from pure numbers to functions depending on the kinematics; 
see also [83]. 

Given these observations, it is plausible that the contribution to the QCD form factor 
remainder with the leading hidden transcendentality should coincide with the correspond- 
ing form factor remainder in MN = 4 SYM theory. Again, the word ‘hidden’ refers to the fact 
that besides harmonic numbers also factors such as TET should formally be attributed 
a degree of transcendentality. A possible starting point to check this conjecture would be 
a comparison between the two-loop dilatation operator of W = 4 SYM theory and QCD, 
see e.g. the recent work [84]. In [53], it was shown that the two-loop dilatation operator in 
the (scalar) SL(2) sector of W = 2 SYM theory is indeed given by the expression in M = 4 
SYM theory plus terms of lower transcendentality, which at least points into this direction. 
It would also be interesting to study the remainder functions of form factors in QCD with 
regard to this point. 
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A Explicit one-loop integral 


In this appendix, we evaluate the one-loop integral (3.7) and show how to obtain the 
explicit result (3.9). 

The integral (3.7), shown here for ? = 1 for simplicity, can be regarded as a certain 
component of the tensor integral 


Dh Viti -diul əl 
(1) 20-08 1 Al bi (uy ee pm İZ umpi +++ ÜZ min) 
I a = — YE ENE A PM IATA, 
25— ni RRO -ny (A) 
AK ye Pi (m "+" Pl pp P2 pega ++ P2 14) 
aa > 1 rr m : 
k=0 


Here, l2 = pı + po — li, sız = (pi + pa)? = 2p) - p2, and the parentheses around the 
Lorentz indices indicate traceless symmetrisation w.r.t. the enclosed indices. The traceless 
symmetrised product built from a single momentum was given in [85] as 


IT(s+n 1) 
7 m= La A ra "o Mponji Hons liissi 27. 
(A.2) 
where $, denotes the group of permutations of n elements and ny, is the metric of D- 
dimensional Minkowski space in the mostly-minus convention. The above expression is 
identical to the ordinary product if the momentum / is on-shell. Using the above result, 
an expression for the numerator of the integral in (A.1) can be found: 
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If the momenta lı and lg are on-shell, the traceless symmetric products on the r.h.s. become 
ordinary products. It is easy to see that in this case the contraction of two of these 
tensor products with rank m + n yields zero unless the two tensor products are related 
by exchanging m and n. This applies to the tensors in the decomposition in (A.1), and it 
allows us to project out the coefficients (1; ya n Of (A.1) by contractions with the tensors 
pya 40000 pik), where k+1=m>+nm: 
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(A.5) 
Here, q = pı + pə with q? = s12 < 0 denotes an off-shell momentum which should not be 
confused with the total momentum used in the main text. These latter expressions follow 
from contractions of (A.3) and (A.2) in which l1, l2 are replaced by pj, pə and £ is replaced 
by q, respectively. The expression for Qm mn, which is manifestly symmetric in m and n 
as expected, was obtained by also making use of the reflection formula for the Gamma 


functions: 
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The numerator in (A.4) is given by contracting (A.3) with pya naD D as pitt) 
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In the integral (A.4), all numerator terms that contain lí or l2 lead to (massless) bubble 
integrals with a single on-shell momentum, which vanish. Hence, the traceless symmetric 
products on the r.h.s. of the above expression can be regarded as ordinary products. More- 
over, lı -l2 can be replaced by pı - p2. These changes do not alter the result of the integral, 
which can thus be expressed as 
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where we also have inserted l2 = pı + p2 — lı and then expanded the obtained binomial 
products. 
The integral that occurs in the above expression is given by 
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as explained below. The diagram denotes the standard scalar bubble integral, i.e. 
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In the case n > 1, the result (A.9) can straightforwardly be obtained by using lı - pı = 
3 (U—(l1 —p1)?) for one of the numerator factors and exploiting the fact that the cancellation 
of l? in the denominator yields a (massless) bubble integral with a single external on-shell 
momentum, which hence vanishes. In the special case n = 0 but m > 1, one shifts the loop 
momentum as lı = İl + p; and expands the power of (hı + pı): p2 in a binomial series, then 
using l1 - pa = (Ê — (li — pə)?) for one of the numerator factors of the resulting integrals. 
The resulting (massless) bubble integrals are either zero since they depend only on one 
external on-shell momentum or they can be evaluated using (see e.g. [86]) 
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(A.11) 
Finally, also the scalar triangle diagram occurs, in particular for n = m = 0 but also in the 


binomial expansions. It can be related to the bubble integral via 
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A resummation of the binomial series yields a result which is in accordance with (A.9) for 
n=0. 


After inserting (A.9) into (A.8), the sums over s and t can be performed, and after 


26 


also applying the reflection formula (A.6), the expression reads 
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(A.13) 
Inserting this result and (A.5) into (A.4), we finally obtain for the coefficients in (A.1): 
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B Remainder densities 


In this appendix, we give the remainder densities up to spin three in the notation introduced 
in (4.28). 
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